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An accurate calculation ol the exciton ground and excited states in AlGaAs and InGaAs coupled 
quantum wells (CQWs) in an external electric field is presented. An efficient and straightfor- 
ward algorithm of solving the Schrodinger equation in real space has been developed and exciton 
binding energies, oscillator strengths, lifetimes, and absorption spectra are calculated for applied 
electric fields up to lOOkV/cm. It is found that in symmetric 8-4-8 nm GaAs/Alo.33Gao.67As CQW 
structure, the ground state of the system switches from direct to indirect exciton at approximately 
5 kV/cm with dramatic changes of its binding energy and oscillator strength while the bright excited 
direct-exciton state remains almost unaffected. It is shown that the excitonic lifetime is dominated 
either by the radiative recombination or by tunneling processes at small/large values of the electric 
field, respectively. The calculated lifetime of the exciton ground state as a function of the bias 
voltage is in a quantitative agreement with low-temperature photoluminescence measurements. We 
have also made freely available a numerical code for calculation of the optical properties of direct 
and indirect excitons in CQWs in an electric field. 
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I. INTRODUCTION 

The electronic and optical properties of quantum well 
structures have been widely investigated in the past 
decade due to their potential applications in electro-optic 
and optoelectronic devices. In recent years there has 
been growing interest in coupled quantum wells (CQWs) 
due to formation of long-lived excitons when these struc- 
tures are placed in an electric field (EF). Intensive stud- 
ies of indirect excitons in CQWs have resulted in their 
electrostatic and optical controlP^ Very recently, CQWs 
have been embedded into Bragg-mirror microcavities and 
there has been found a special type of voltage-tuned exci- 
ton polaritons which can be used for optical nonlinearities 
and polariton lasing achieved at much lower threshold 
powers P 

A CQW structure consists of two quantum wells sep- 
arated by a barrier layer. For a sufficiently thin barrier, 
the tunneling of carriers through the barrier makes the 
two wells electronically coupled to each other. As a re- 
sult, an electron (hole) can either reside in one of the two 
wells, or its wave function (WF) is distributed between 
both wells. In the case of Coulomb bound electron and 
hole residing in the same well, they form a direct exciton. 
If however they are located in different wells, an indirect 
exciton is created. 

In a symmetric CQW structure with no EF applied, 
formerly degenerate single-particle states split, owing to 
the tunneling through the middle barrier, into doublets 
with symmetric and antisymmetric states in each. Since 
only transitions between states having the same parity 
are optically allowed, the Coulomb-coupled electron-hole 
(e-h) pairs form excitonic states which are optically ei- 
ther bright or dark. An EF, being applied in the growth 
direction, breaks down the symmetry of the system mak- 



ing all these excitons bright. In fact, single-particle states 
experience with EF a transition from states with well de- 
fined parity to the ones with the electron (hole) located 
in one of the two wells, thus forming direct and indirect 
combinations of uncorrelated e-h pair states. These dif- 
ferent pair states are Coulomb coupled with each other 
and form an exciton in which direct or indirect pair 
can dominate. In particular, with increasing EF, the 
exciton ground state (GS) undergoes a transition from 
bright direct exciton to indirect exciton which has a much 
weaker optical activity. The exciton radiative lifetime in- 
creases due to a reduction in the spatial overlap between 
the electron and hole WFsPHS The exciton binding en- 
ergy, in turn, reduces owing to an increased e-h sepa- 
ration. It has also been found that the electronic cou- 
pling between quantum wells considerably enhances the 
quantum-confined Stark effect in CQW structuresf^EZl 
The tunneling effect is also enhanced with the E F allow- 
ing the carriers to leak out of the system ESU2I This can 
lead to a considerable shortening of the photolumines- 
cence decay timeP3 All these properties of CQWs make 
them a much richer system compared to single quantum 
wells. 

Excitonic states in CQWs in the presence of EF have 
been intensively studied in recent years. Different the- 
oretical approaches have been used ranging from varia- 
tional method a 21 " 2 ' to direct diagonalizations in which 
the exciton WF is expanded into a large basis^ or the 
Schrodinger equation is discretized in the momentum 
spaced In this paper, we present a more accurate and 
straightforward way for solving the Schrodinger equation 
for an exciton in a CQW structure. Expanding the exci- 
ton WF into e-h pair states we solve in the real space a 
system of differential equations for the exciton in-plane 
motion, using the shooting method, here generalized to 
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a matrix form. The e-h pair basis states are calculated 
exactly using the analytical form of the electron (hole) 
WF in a uniform EF. We present the full calculation of 
exciton bound and (discretized) continuum states as well 
as the absorption spectrum of a CQW in an applied EF. 
We study the EF effect on the exciton binding energy 
and lifetime in 8-4-8 nm GaAs/Al .33Gao.67As symmetric 
CQW structure, both for the ground and excited states, 
and demonstrate a direct-to-indirect crossover of the ex- 
citon ground state with increasing EF. An example of an 
asymmetric 10-4-10 nm InGaAs CQW having different In 
content in the left and right QWs is also given, demon- 
strating our calculation of the electron and hole energies 
of quantization and the exciton transition energies and 
oscillator strengths as functions of the applied EF. 



II. FORMALISM AND NUMERICAL METHOD 

Let us consider a symmetric GaAs/Al^Gai-^As CQW 
which consists of two GaAs QW layers separated by a 
thin Alj.Gai_3.As barrier and surrounded on both sides 
by thick barriers of the same kind. In this paper we 
mainly concentrate on a CQW structure which has been 
intensively used in a series of experiments,'^^' tak- 
ing the barrier and well widths to be Lf, = 4nm and 
L w = 8 nm, respectively, and the barrier concentration 
of Al to be x = 0.33. The electric field F is applied in 
the growth direction. For a different CQW in presence 
of an EF simulations can be performed using our online 
available numerical code^ (see Sec. HID below for more 
details). We are interested in optically allowed transi- 
tions in such a system and thus consider excitonic states 
with zero in-plane and angular momenta only. In the 
effective mass approximation, the excitonic Hamiltonian 
can be divided into three parts: the first two, H e and Hh, 
take into account the electron and hole quantization in 
heterostructure potentials V e and Vh and the third one, 
Hx, is responsible for the electron- hole (e-h) in-plane 
relative motion and Coulomb binding: 

H{z e ,z h ,p) = H e (z e ) + H h (zh) + H x {z e , z h , p) + E g (1) 

with 
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where z e (^h) is the electron (hole) coordinate in the growth 
direction, p the coordinate of the e-h relative motion 
in the QW plane, the background dielectric constant 
which we assume to be z-independent, m e the electron 
effective mass, and E g the bandgap of the well mate- 
rial (GaAs). Owing to the strong QW confinement, the 
heavy-hole subband is split off considerably and can be 
approximated by an anisotropic effective mass using the 
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where mh is the hole effective mass in the growth direc- 
tion, p the exciton in-plane reduced mass, and mo the 
free electron mass. We assume a rectangular form of the 
heterostructure confinement potentials, 



V e ,h{z) 
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inside the wells, 
outside, 
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and similar step-like profiles for the electron and hole 
effective masses in the growth direction. 

In our calculation, we have used the following param- 
eters: the background dielectric constant 24 Eh = 12.5, 
the energy-band offset ratio 2 ^! y e : V h = 65 : 35. The 
band gap discontinuity at the GaAs/Al^Gai-^As inter- 
face is linearly approximated^ as 1.247-xeV. The Kohn- 
Luttinger parameters for pure GaAs and AlAs are ob- 
tained from Ref.1321 m e = 0.0665 mo, 71 = 6.79 and 
72 = 1.92 in GaAs (giving the hole mass ran = 0.34 mo); 
m e = 0.15 m , 71 = 3.79 and 72 = 1.23 in AlAs (giv- 
ing m/j = 0.75 mo). The AL-Gai-^As alloy parameters 
were linearly interpolated between those of GaAs and 
AlAs. In particular, for the content x = 0.33, we have 
used m e = 0.094 m and m h = 0.48 m in Alo.33Gao.67As 
layers. The in-plane reduced mass in the barrier lay- 
ers (0.057 mo) is different from that in the well layers 
(0.042 mo). However, because of a very small probability 
for the exciton ground state to find the carriers in the 
barrier, we take for the in-plane reduced mass the GaAs 
value of p — 0.042 mo. 



A. Electron and hole single-particle states 

To solve the excitonic Schrodinger equation with the 
full Hamiltonian Eq. (JlJ we first prepare a basis of 
single-particle states which satisfy the following one- 
dimensional equations: 



(7) 



To do so, we use the advantage of the analytic form of 
the electron and hole WFs in the rectangular confinement 
potentials Eq. Q and uniform electric field F. In each 
layer of the CQW structure, the electron WF is given by 
a superposition of two Airy functions^ 



r(z) = a fc Ai(0 + & fc Bi(£) , 
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and the index k labels the heterostructure layers (from 
left to right) taking integer values from 1 to 5. The elec- 
tron eigenenergy E e and five pairs of coefficients (a,k,bk) 
in Eq. Q are found from four pairs of boundary condi- 
tions (BCs) on heterostructure interfaces and two BCs at 
z — > ±00. 

The interface BCs following from Eqs. (pi) and ([7| are 
the continuity of ip e (z) an d m~ 1 (z)di/j e (z)/oz. The other 
two BCs take into account the possibility for the electron 
to tunnel through the barrier and escape from the system 
to the side of the CQW structure where the applied EF 
gradually lowers the potential. In that area, the solution 
is given by a wave propagating away from the system. For 
the electron WF and F > 0, this outgoing BC at z — > —00 
yields b\ = —%a\ which follows from the specific combi- 
nation of the Airy f unctions Ai(£) — «Bi(£) producing an 
outgoing waveP^ESl At the same time, the electron cannot 
escape to the other side of the structure where the poten- 
tial gradually increases, and thus the other BC, typical 
for bound/localized states, is ip e (z — > +00) = 0, giving 
65 = due to the asymptotics of the Airy functions.^ 
The form of the WF and the BCs for the hole are found 
in a similar way, taking into account that the potential 
grows in the opposite direction. The secular equation fol- 
lowing from all ten BCs determines discrete eigenvalues 
of Eq.Q, 
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which are the complex energies of electron/hole resonant 
states, also known in the literature as Siegert states.^ 
The real part of the eigenvalue, E^ h ' , is the energy po- 
sition of the electron (hole) j-th resonant level, while the 
imaginary part gives its tunneling linewidth. 

The WF of any resonant state having a finite linewidth 
is essentially complex, i. e. it cannot be made real by any 
uniform phase shift. Also, its amplitude grows exponen- 
tially to the outside area to which the partic le ca n es- 
cape and thus has to be normalized to its flux! 35 * 36 ^ This 
normalization includes a divergent volume integral and 
a compensating surface term. For the values of the EF 
considered in this paper, the calculated linewidths of the 
electron and hole states of interest are always small com- 
pared to their energies of quantization. Similarly, the 
imaginary parts of the WFs are small compared to the 
real ones and can be dropped. The normalization condi- 
tion is then taken in a form 



$f(z)] 2 dz 
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where = Re(V>^'' t ). 



Eq.(ll| 



The limits of integration in 
and Zmax, taken to be the same for elec- 
tron and hole, are two distant points on both sides of 
the CQW where the WFs decay considerably before they 
start to grow exponentially owing to the carrier tunnel- 
ing, so that the surface terms are minimized and can be 
dropped leaving in the normalization only finite-volume 
integrals. 
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FIG. 1: Energies of two electron states (a) and four hole states 
(b) in 8-4-8 nm GaAs/Alo.33Gao.67As CQW as functions of 
an applied electric field. Insets zoom in particular spectral 
regions with anticrossing. 



Different electron and hole subbands (Ef for i = 1,2 
and Ej for j = 1-4) calculated in the presence of the 
EF are shown in Fig.Jl] The corresponding WFs, t/>f 
and ipj, are illustrated in Fig.[2ja)-(d) for a few dif- 
ferent values of the EF. At zero field, the GS and the 
first ES have, respectively, symmetric and antisymmet- 
ric WFs, see Fig.llja). With increasing EF, the WFs 
become asymmetric, and the WF maxima for the elec- 
tron and hole GSs move in opposite directions. The GS 
and the first ES for the same carrier are also confined 
in different QWs. This happens to both carriers already 
at F = 2kV/cm, see Fig.[5Jb). As the EF grows fur- 
ther, the GS-ES splittings increase almost linearly with 
F (Fig. [I]), and the hole ES jumps from the left to the 
right QW [Fig.ptc)]. This corresponds to an anticrossing 
of hole ES subbands which takes place at 27.8kV/cm, see 
the inset in Fig.JTJb). Such an anti crossin g behavior was 
also found in previous calculations. 1 ^!! The same hap- 
pens to the electron ES at a much higher EF. Figure [2jd) 
shows a modified ES WF for F = lOOkV/cm which is a 
precursor of a similar transition for the electron. 
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FIG. 2: (a)-(d) Wave functions and energy levels of electron 
and hole ground and excited states for different values of the 
electric field F. Grey lines show CQW heterostructure poten- 
tials, using Eg = 1.519 eV for the band gap. (e),(f) Oscillating 
tails in the electron and hole wave functions. 



Far from the CQW structure, we observe tiny oscil- 
lations in the WF of the electron (hole) GS in the left 
(right) barrier, see Fig.[2je,f). This oscillatory behav- 
ior, typical for freely propagating particles, occurs due to 
the lowering of the potential by the EF, so that an elec- 
tron (hole) can escape to the left (right) barrier. Since 
nihVh > rrieVet the frequency of the oscillations for the 
hole is larger than that for the electron. At the same 
time, the amplitude of oscillations is also larger for the 
hole since the hole tunneling is stronger due to Vh <V e . 
Such oscillations become more dramatic and start ear- 
lier (closer to the CQW center) for some of higher ESs. 
For example, for F = lOOkV/cm the hole state h2> ex- 
ibits huge oscillations clearly seen in Fig.[2]jd). However, 
the next hole ES hA which has the dominant contribu- 
tion to the direct exciton state discussed below has only 
tiny oscillations, similar to those in Fig.^f), and a small 
tunneling rate, comparable to that of the hole GS. 



B. Excitonic states: Multi-sublevel approach 

We calculate excitonic states in a CQW structure in 
the presence of an EF, expanding the exciton WF into a 
finite set of e-h pair states: 



$(z e , Z h ,p) = ^ ®n(Ze, Zh)<f>n(p) 



(12) 
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where 



$ n {z e ,z h ) = i>i(ze)ipj(.Zh), n = {i,j), (13) 



and (z) are the electron and hole wave functions cal- 
culated in the presence of EF and heterostructure po- 
tentials, see Sec. |II A| The Schrodinger equation for the 
exciton, — Ex^, then takes the form 
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and the energies of pair states given by 
E (o) 
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In our calculation of the exciton states in 8-4-8 nm 
GaAs/Alo.aaGao.erAs CQW for EFs up to F = 25 kV/cm 



it was sufficient to restrict the basis in Eq. (12 1 to four 
e-h pair states (N — 4), keeping only the GS and the 
first ES for electron and hole. We label these four basis 
states as elhl (n = 1), elh2 (n = 2), e2hl (n = 3), and 
e2h2 (n — 4), where e(h) stands for an electron (hole) 
and the numbers 1 and 2 refer to the single-particle GS 
and first ES, respectively, while n is the unified pair in- 
dex introduced in Eq. ( fl3| . Other QW structures may re- 
quire higher ESs to be taken into account. These are also 
needed in our case when a stronger EF is considered. In 
particular, higher quantized levels for the hole h3 and M 
[see Fig. 1(b)] are taken into account for F > 25kV/cm 
and F > 70kV/cm, respectively, and these states have a 
major contribution to the direct exciton WF. 

We calculate the exciton transition energy Ex and 
the in-plane components of the WF, (f> n (p), by solving 
the matrix differential equation ( 14 ) numerically. To do 



so, we introduce a matrix generalization of the shoot- 
ing method applying the latter to a system of coupled 
differential equations. The shooting method transforms 
a boundary-value problem like Schrodinger's equation 
with BCs to an initial- value problem in which one of the 
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FIG. 3: Matrix elements V nm of the Coulomb potential cal- 
culated between different e-h pair states at _F=4kV/cm. 



boundary values (in the present case the WF at p — > oo) 
is taken as a starting point. The boundary value(s) on 
the other side (at p = 0) is then used to find the eigenen- 
ergies. The BCs follow straightforwardly from Eq. (14) 



and the asymptotics of the Coulomb matrix elements 



V nm (p). At large distances V nm (p) 



7(M> 



while at small distances the potentials V nm (p) have log- 
arithmic dependence, as is clear from Fig. 3. Therefore 
for bound states 



4> n (p ->• oo) = A n p Sn e~ 
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Discretizing Eq. ( 14 ) on a finite grid, a numerical solu- 
tion in the area < p < R is generated iteratively using 
a finite difference scheme. In particular, a second-order 
scheme which we have used in our calculation brings 



Eq. ( 14 1 to the form 



4> n (p-Ap) = -<f> n (p + Ap) 



2p + Ap 
2p- Ap 
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E F nm{p)<P m {p) , 



(20) 

where Ap is the discretization step and the matrix 
Fnm(p) depends on the Coulomb interaction and on 
a finite-difference representation of the kinetic term 
Eq. ( p~5| ) . The WF amplitudes A n in the starting values 
(f) n (R) in Eq. (18) are the unknowns to be found along 
with the eigenvalue Ex- Going to the very last point 
p = and using the boundary condition Eq. ( 19 1 pro- 



duces a homogeneous matrix equation of the form 



M nm {E x )A 

r, 



(21) 



in which M nm (Ex) depends solely on the exciton energy 
Ex (and not on p any more), and thus the energy eigen- 
values are determined by 



det|M nm (£ x )| =0. 



(22) 



For small values of the EF, the electron and hole GS-ES 
splittings are smaller than the exciton Coulomb energy 
(compare Figs. 1 and 3) and thus several bound states 
[having the asymptotics given by Eq. ( 18 )] can always 
be found in the system. However, as the EF grows, the 
Coulomb energy of the exciton ESs is getting smaller than 
the e-h pair splitting energies and thus some of these ex- 
citon states become unbound. Since the unbound states 
have energies Ex > E^ , at least for some of their radial 



components the asymptotics Eq. ( 18 ) is no longer valid 



and a proper treatment of the excitonic continuum is re- 
quired. This task is outside the scope of the present paper 
which mainly concentrates on exciton bound states. Nev- 
ertheless, some effects of the continuum and in particular 
its influence on the excitonic absorption spectrum can be 
taken into account, in a first attempt, by restricting the 
exciton in-plane motion to a large circle of radius R and 
in this way discretizing the continuum. The asymptotic 
boundary conditions Eq. (18) are now replaced by 



(f> n (R) = 0, 4> n {R - Ap) = A, 



(23) 



where the new amplitudes A n satisfy the same Eq. (21) 
with matrix M nm (E) being redefined accordingly. 

It is convenient to normalize the radial components of 
the WFs introducing expansion coefficients C n : 



4>n(p) = C n (f> n (p) , 



(24) 



where 4> n (p) is normalized to 2-7T \(j) n \ 2 pdp = 1, and 
therefore 



(25) 



due to orthogonality of the e-h pair states (fr n (z e , z^) and 
normalization of the total exciton WF ^>(z e , Zh, p)- 

Finally, for each excitonic state, the oscillator strength 
per unit area is calculated as^ 



2m E x \d c 
H 2 



ty(z,z,p = 0) dz 
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where d cv is the basic dipole matrix element between the 
valence and conduction bands, and the overlap integral 
in Eq. ( 26 1 accounts for the spatial distribution of the 



excitonic recombination. The exciton radiative linewidth 
is then given by 



e b m c 
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(27) 



where c is the speed of light. 



III. RESULTS AND DISCUSSION 

In the presence of EF, the energy levels of the elec- 
tron and hole single-particle states in a CQW structure 



experience Stark shifts (F) w E{^' (0) T dF/2 as 

demonstrated in Fig. 1, where d — L w + Lf, is the center- 
to-center distance between the QWs. Concentrating on 
these two lowest levels for the electron and two for the 
hole we are thus dealing with four e-h pair states. For tw o 
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of them, elh2 and e2M, the energies E^l [see Eq. ( IT )] 
remain almost unaffected by the EF as the Stark shift 
for the electron is compensated by that for the hole, 
while the other two pair states, elhl and e2h2, have 
twice larger Stark shifts than the single-particle states: 



E[°l(F) w E\"i(0) T dF. Even though our CQW is sym- 
metric, at nonzero EF the ground and the first excited 
single-particle states are localized in different wells of the 
CQW structure, see Fig. 2(b). That is why the pair states 
with electron and hole in the same QW (n = 2 and n = 3) 
are electrically neutral and can be called direct states 
while the other two, n — 1 and n = 4, with electron and 
hole in different QWs have nonzero dipole moment and 
are called indirect states. As for excitonic states, they 
are, strictly speaking, neither direct nor indirect, since 
the exciton WF is always a combination of different pair 
states. 
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A. Crossover from direct to indirect exciton 

In the absence of the Coulomb interaction, all four pair 
states are nearly degenerate at zero EF: Due to a rather 
weak tunneling through the middle barrier of the CQW, 
the energy splittings are small compared to the Coulomb 
energy (2.3 meV and 0.05 meV splitting for electron and 
hole, respectively, versus 4-9 meV binding energy). The 
Coulomb interaction splits further and strongly mixes the 
direct and indirect e-h pair states. To understand these 
splitting and mixing, let us consider the Coulomb matrix 
elements in more detail. 

Figure 3 shows an example of the Coulomb matrix el- 
ements for F — 4kV/cm, though this picture does not 
change much when the EF increases or decreases. The 
direct and indirect pair states have different charge sepa- 
ration, and thus their diagonal Coulomb matrix elements 
are also quite different: Potentials V22 ~ V33 for the di- 
rect pairs are a few times stronger than V\\ ~ V44 for 
the indirect ones. This difference brings in a consider- 
able splitting between the direct and indirect pair states 
at F — 0: The indirect doublet is found almost 5meV 
above the direct one, see Fig. 4. Among off-diagonal el- 
ements, the largest are Vi 3 ps V24, due to a considerable 
overlap between the electron ground and excited states. 
All other matrix elements are two-three orders of mag- 
nitude smaller, because of a much larger effective mass 
of the hole and consequently much smaller overlap inte- 
grals. All off-diagonal elements drop quickly at large p 
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FIG. 4: Electric field dependence of the optical tran- 
sition energy Ex for different exciton states in 8-4-8 nm 
GaAs/Alo.33Gao.67As symmetric CQW structure, calculated 
using the exciton confinement radius R = 200 nm (a) and 
800 nm (b). The circle area is proportional to the exciton 
oscillator strength /. 



due to the orthogonality of WFs. 

The Coulomb coupling matrix elements V13 ~ V24 are 
responsible for the mixing of direct and indirect pair 
states and in particular for a crossover of the ground ex- 
citon state from direct to indirect type as the EF grows. 
Owing to this off-diagonal coupling, the Stark red-shifted 
indirect state elhl (n = 1) has a remarkable anticrossing 
with the direct pair state e2hl (n — 3), weakly depen- 
dent on the EF. This anticrossing takes place at about 
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FIG. 6: Radial components 2-k p\4> n (p)\ 2 of the exciton wave 
function calculated for the ground state X-GS (a-d) and ex- 
cited state X-ES (e-h), for different values of the electric field. 



FIG. 5: Probability distributions J Q °° |$(2: e , z^, p)\ 2 2npdp cal- 
culated for the exciton ground state X-GS (a-e) and excited 
state X-ES (f-j), for different values of the electric field. 



F = 5kV/cm and is clearly seen in Fig.|4j Indeed, the 
exciton ground state (X-GS) experiences a crossover from 
direct to indirect exciton: The oscillator strength of the 
X-GS has its maximum at F — and then drops quickly 
with increasing EF as seen in Fig. [4] Since all other ma- 
trix elements including V12 are generally much smaller, 
the other direct pair state elh2 (n = 2) remains unaf- 
fected and is only Coulomb shifted by the diagonal ele- 
ment V22- Although it is strongly coupled via V24 to the 
other indirect state e2h2 (n = 4), the latter is Stark blue- 
shifted and thus significantly detuned from elh2 having 
a minor effect on it. As a result, the energy position of 
the brightest exciton excited state (X-ES), which has the 
maximum oscillator strength in the excitonic spectrum, 
remains practically unchanged. 



The excitonic states shown in Fig. 4(a) and (b) are cal- 
culated with different in-plane exciton confinement ra- 
dius, R = 200 nm and 800 nm, respectively [see Eq. (23)] 



In the latter case the excitonic continuum has a much 
finer discretization that makes more clear which states 
belong to the continuum and which are the true bound 
states having more or less isolated positions, weakly de- 
pendent on R. For example, the 2S* and 3S states of 
the indirect exciton are clearly identified in Fig. 4(b): 
They lie just below the discretized continuum onset and 
are down-shifted with F almost parallel to the exciton 
ground state. Higher excited states of the direct exciton 
are also well seen in Fig. 4: They are deep in the contin- 
uum (7-8 meV above the brightest direct state) and are 
weakly dependent on the EF. 

Let us consider the properties of the X-GS and X-ES, 
and in particular the direct-to-indirect (D-I) crossover in 
more detail. The contour plots in Fig. [5] show localization 
of the X-GS and X-ES across the CQW structure, for dif- 
ferent values of F. At zero EF, owing to the symmetry of 
the system, both states have two identical maxima on the 
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FIG. 7: Coefficients C n of the expansion of the exciton wave 
function into e-h pair states, calculated for the ground state 
X-GS (a) and excited state X-ES (b) as functions of the elec- 
tric field F. 



main diagonal z e — Zh [Fig.[5ja) and (f)], which refer to 
the direct nature of both excitonic states. With increas- 
ing EF one of the two peaks becomes smaller and then 
vanishes, see Fig.[5]jb,c) and (g,h). The states become 
asymmetric having both carriers localized either in the 
right QW (in X-GS) or in the left QW (in X-ES). Further 
increase of the EF up to F = 6 kV/cm leads to the X-GS 
switching from direct to indirect state: The peak moves 
away from the main diagonal towards the bottom right 
corner, see Fig.gc)-(e). This result is in good agreement 
with previous theoretical findings^EH anc j experimental 
observations.^ At the same time the X-ES remains un- 
changed. The dominant hole component in the X-ES 
changes from h2 to h3 and from /i3 to bA when the sys- 
tem passes through the anticrossings of the hole levels 
at F = 27.8kV/cm and 74kV/cm, respectively. The X- 
ES remains the brightest state in the excitonic spectrum 
having at the same time a very weak dependence on the 
EF that emphasizes its direct nature. 

The D-I crossover of the X-GS is demonstrated also 
in Fig. 6 where different radial components 4> n (p) of the 
X-GS and X-ES WFs are plotted. At F = 2kV/cm the 
direct e-h pair state n = 3 has the dominant contribution 
to the X-GS. This state is strongly coupled to the indi- 
rect pair state n — 1 via the Coulomb matrix element 
Vi3, as discussed above. As a result of this anticrossing, 
(j>i grows and 03 reduces with the EF. Nothing similar 
happens to the X-ES. The latter is always dominated by 
the n = 2 pair state which is coupled to the significantly 
detuned n — 4 state. Therefore, only a minor contribu- 
tion of the (f>4 to the X-ES can be seen in Fig. 6(f)-(h). 
The picture is completely different for smaller EFs. At 
F = the symmetric n = 1 pair state is Coulomb coupled 
to the other symmetric state (n — 4) producing bright 
X-GS and dark X-ES. These symmetric pair states do 
not interact with the antisymmetric states (n = 2 and 3) 
which are, in turn, coupled to each other. The transition 
from such symmetric coupling to the above considered 
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rate r_R of the exciton ground state X-GS (solid lines) and 
excited state X-ES (full squares) as functions of the electric 
field F. Dashed lines are the single state basis calculation of 
the X-GS. 



D-I coupling has a rather narrow interval of small val- 
ues of the EF and involves interaction of all four e-h pair 
states. In fact, for F — 0.1 kV/cm one can see in Fig. 6(a) 
that all four components have comparable contributions 
to the X-GS and X-ES WFs. Figure 7 summarizes our 
analysis showing the field dependence of the e-h pair am- 
plitudes C n introduced in Eq. ( 24 ) . It demonstrates the 



prominent D-I crossover in the X-GS, a much weaker D-I 
coupling in the X-ES, and a very quick transition from 
symmetric coupling to D-I coupling, seen in the WFs of 
both X-GS and X-ES. 

Figure 8 shows the field dependence of the optical 
transition energies Ex and binding energies E), of X-GS 
and X-ES, as well as their in-plane exciton Bohr radii 
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r B = VV) and the radiative linewidths Tr. The l atter 
are calculated via Eq. (27 1 assuming d cv = 0.6 nmP^22l 



While all these parameters for the X-GS change dramat- 
ically when the EF increases from 2kV/cm to 8kV/cm, 
the X-ES remains practically unaffected. The X-GS en- 
ergy exhibits a considerable Stark shift, for F > 6kV/cm 
almost linear in F [Fig.[8ja)], due to the electron- hole sep- 
aration typical for the indirect exciton. The X-ES in turn 
has a very weak field dependence due to a much smaller 
dipole moment of the direct exciton, but in a larger range 
of the EF values, the X-ES transition energy is also red- 
shifted (by 12meV at F = lOOkV/cm), as a result of 
an EF-induced electron-hole separation within the same 
QW. 

The X-GS binding energy, E b = — E x , drops from 
8meV down to 4meV [Fig.|8jb)], as a result of the tran- 
sition from direct to indirect Coulomb coupling. For the 
X-ES the binding energy is defined as the energy distance 
from X-ES to its own continuum onset: E\, = E^ — Ex- 
This is done because the lowest-energy pair state, el hi, 
has negligible (for F > lkV/cm) contribution to the 
X-ES, and thus this exciton state remains bound even 
though its energy Ex is higher than the first contin- 
uum onset e[ 0) . The X-ES bindin g energy is almost un- 
changed, as the EF has no effect on the direct nature of 
this state. The Bohr radius [Fig.|8jc)] is fully correlated 
with the binding energy, increasing with the EF almost 
by a factor of 2 for the X-GS and showing no change for 
the X-ES. 

The D-I crossover is accompanied by a dramatic de- 
crease of the radiative linewidth of the X-GS (propor- 
tional to its oscillator strength), see Fig.lstd). This hap- 
pens because the EF makes the electron-hole separation 
larger reducing the overlap of the electron and hole WFs. 
The X-ES, in turn, being strictly dark at F — 0, becomes 
bright in a finite EF, and its linewidth quickly increases 
with the EF up to the half of the maximum linewidth of 
the X-GS. Further increase of the EF does not change 
the X-ES radiative rate much but Tr experiences some 
influence of higher ESs, see fluctuations in Fig.[8^d). 

To reproduce some previous simulations^ and to com- 
pare them with our full calculation, we have restricted 
our basis to the electron and hole GSs only, neglect- 
ing any ESs quantized in the growth direction. In other 
words, we have taken into account only one e-h pair state 
elhl, choosing N — 1 in the expansion Eq. (12 1. This 



single-state basis (SSB) calculation shows considerably 
different results compared to the full calculation in all 
four plots in Fig.[8j see dashed curves. The reason for 
such difference is obvious: Figure [7][a) clearly demon- 
strates the importance of taking into account higher e-h 
pair states for proper description of the D-I crossover of 
the X-GS, and in particular the role in such crossover of 
state elhl which is missing in the SSB calculation. This 
state is a direct pair state (at least for F > 0.5kV/cm), 
so that neglecting it in the X-GS calculation, as done 
e.g. in Ref. 1211 underestimates the exciton binding en- 
ergy by a factor of 1.5 and the exciton radiative linewidth 
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FIG. 9: (a) Exciton energies and oscillator strengths (cir- 
cle area) and (b,c) electron and hole subbands in asymmet- 
ric Ino.08Gao.92As/GaAs/Ino.1Gao.9As (10nm/4nm/10nm) 
CQW structure as functions of the applied electric field F. 
The inset shows CQW band structure and direction of the 
electric field, as well as the electron and hole wave functions 
and energy levels of the ground and first excited quantized 
states, for F = 12kV/cm. 



by an order of magnitude. Nevertheless, at larger EFs 
(F > lOkV/cm) the SSB model adequately describes 
the properties of the X-GS as such an indirect exciton 
state is strongly dominated by elhl. 

Finally, Fig. 9 shows our simulation for an asymmetric 
Ino.08Gao.92As/GaAs/Ino.1Gao.9As (10nm/4nm/10nm) 
CQW structure used in Ref. [HI As expected, the en- 
ergy spectra are asymmetric with respect to the EF di- 
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can have some effect on the exciton tunneling, its life- 
time is strongly dominated by the radiative channel as 
can be seen in Fig.[l0j In particular, the tunneling life- 
time t t = h/(Tl + Vf) is much longer than the radiative 
one up to F = 80 kV/cm. The D-I crossover of the X-GS 
is accompanied by a monotonous growth of its radiative 
lifetime, in accordance with Fig. 8(d). Indeed, a direct ex- 
citon has a short lifetime because the carriers are in the 
same well, so that they can easily recombine. Increasing 
the electron-hole separation leads to a dramatic increase 
of the radiative lifetime. The probability of the elec- 
tron and hole tunneling also increases with EF. At some 
point the tunneling time becomes comparable to the ra- 
diative lifetime and then starts to dominate. We have 
also compared the calculated radiative lifetime for the 
X-GS with the experimental results taken from Ref. [501 
where the low-temperature excitonic photoluminescence 
was measured in 8-4-8 nm GaAs/Alo.33Gao.67As CQWs. 
The inset to Fig.JTUJ demonstrates a quantitative agree- 
ment between the experiment and the present theory. 



C. Absorption spectrum 



rection. In the right part of the exciton energy spectrum 
[Fig. 9(a)], at around 16kV/cm, a Coulomb-induced anti- 
crossing is observed, which is similar to that seen in Fig. 4 
at F = 5kV/cm. The physical mechanism which causes 
this anticrossing is essentially the same as in the case 
of symmetric CQWs, but the anticrossing takes place at 
much higher values of the electric field. This is because at 
around F = 12kV/cm (4-5kV/cm below the anticross- 
ing) the applied EF almost compensates the asymmetry 
in the conduction band structure, and thus at this value 
of F the properties of the asymmetric CQW can resem- 
ble those of the symmetric CQW at F = 0. In fact, 
a repulsion of the electron subbands is seen in Fig. 9(b) 
as well as a formation of symmetric and antisymmetric 
electronic states, as is clear from the inset in Fig. 9(a). 



B. Indirect exciton lifetime 

A CQW exciton can escape from the system using the 
following two major channels: It can either recombine by 
emitting a photon or the electron and/or hole can tunnel 
through the external barrier with the help of the EF. 
We concentrate here on the X-GS only, and combining 
both channels together, the total exciton lifetime takes 
the form 



1 



1 



(28) 



where r R = K/{2Tr) is the exciton radiative lifetime. As 
for the tunneling time r r , we take into account the low- 
est pair state only. This is a valid approximation because 
when other pair states contribute to the X-GS and thus 



We also calculate the exciton absorption coefficient in 
a 8-4-8 nm GaAs/Al .33Ga .67As CQW at different fre- 
quencies of the incoming li ght. Using the Lorentzian 
model of absorbing oscillator j 37 l 4 ° l and leaving out a com- 
mon prefactor, the absorption of the exciton with zero 
in-plane momentum takes the form 



«M = E r 



R,v 



R,u 



V 2 

1 R,v 



(29) 



where the index v labels all possible excitonic states cal- 
culated in the theory, and E v and stand for their 
energies and radiative linewidths. 

Figure fllTa) shows the absorption spectrum calculated 
using Eq. ( |29[ ) and in-plane exciton confinement radius 
R = 800 nm. All lines in the absorption have very nar- 
row radiative widths (< 0.1 meV) which are calculated 
according to Eqs.(26) and ([27]). The Lorentzian model 
however has an obvious artefact: Although the spec- 
trum properly reproduces the linewidths, all lines have 
the same peak height, and the fact that not all of them 
are seen in Fig. 11 'a) is only due to the resolution of the 
plot. To improve on this and also to take into account 
the effect of inhomogeneous line broadening in realistic 
CQW structures, we make a Gauss convolution of the 
spectrum: A(lu) = J_ ce(u> )g(oj — u )duj with a nor- 
malized Gauss function g(u) = (Ay/n)~ 1 e~ u2 / A ~ . The 
convoluted spectrum with the full widt h a t half maxi- 



b). All lines 
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mum 2v In 2A = 0.2 meV is shown in Fig. 
now have almost the same width but their peak maxima 
now reflect the optical strength of the corresponding ex- 
citon states. 

The two lowest excitonic states, X-GS and X-ES, are 
well resolved in the spectrum up to F — 5kV/cm. Then 
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FIG. 11: (a) Electric field dependence of the full excitonic 
absorption spectrum, (b) The same spectra convoluted with 
a Gaussian function with 0.2 meV full width at half maximum. 



the indirect exciton X-GS loses its optical activity. The 
bright direct X-ES line superimposes with higher ESs 
and discretized continuum of the indirect exciton, all line 
merging up together at higher EFs. 



input, the user has an option either to choose a symmetric 
AlGaAs CQW inserting required structural parameters 
(Lf,, L w , and x), or to take other semiconductor CQW 
by inserting its structural and material parameters. The 
value of the EF is also needed for the input. In this nu- 
merical code, we use four basis states for the exciton WF, 
taking into account only the ground and the first excited 
states for the electron and hole. The exciton confinement 
radius is set to R — 200 nm. The transition energies and 
the oscillator strengths of exciton excited states are also 
included in the output file. The absorption spectrum is 
calculated with 0.1 meV FWHM Gaussian convolution. 



IV. CONCLUSIONS 

We have studied the effect the electric field has on ex- 
citonic states in AlGaAs and InGaAs CQWs. To do this 
we have developed an efficient numerical approach which 
is based on expanding the exciton wave function into un- 
corrclated electron-hole pair states and solving in the real 
space a matrix Schrddinger equation for the CQW exci- 
ton. Using this approach we have calculated the energies 
and the wave functions of exciton states in the CQW 
and studied their optical properties in the presence of 
electric field, by addressing such important parameters 
of the exciton as its binding energy, Bohr radius, radia- 
tive and tunneling times, and optical absorption spec- 
trum. While we are able to calculate a large number of 
exciton states, we have concentrated on two most im- 
portant ones, the exciton ground state X-GS and the 
brightest excited state X-ES. We have shown that the 
Coulomb coupling between direct and indirect pair states 
leads to a prominent effect in the presence of the electric 
field: a direct-to-indirect crossover of the X-GS. At the 
same time, the properties of the X-ES remain almost un- 
changed. Finally, we have calculated the exciton lifetime 
which consists of two main components, radiative and 
tunneling times, and shown that the latter reduces the 
total lifetime at higher electric fields. 



D. Numerical code 

We have provided a freely available on-line software^ 
which calculates the optical properties of the direct and 
indirect excitons in a CQW structure in a perpendicular 
EF. The software has a user-friendly interface and pro- 
duces an on-screen output of the requested calculation as 
well as a pdf version of the same data and plots. For the 
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